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We present a quantum critical behavior of the renormahzed single-particle Wannier function, 
calculated in the Gutzwiller correlated state near the insulator-metal transition (IMT) for cubic 
lattices. The wave function size and its maximum, as well as the system energy scale with increasing 
lattice parameter R as R". Such scaling is interpreted as the evidence of a dominant role of the 
Coulomb repulsion. Relation of the insulator-metal transition lattice-parameter value R = Rc to 
the original Mott criterion is obtained. The method is tested by comparing our results with the 
• exact approach for the Hubbard chain. 
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Insulator-metal transition (IMT) of the Mott-Hubbard type represents one of the central problems in quantum 
I— I, physics, since it exemplifies a transition from well defined atomic states to delocalized (Bloch- or Fermi-liquid- type) 
I states in a solidi and in other systemsi^ The model system is a lattice of spin-1/2 fermions with one particle per 
I [ atomic site. Mott expressed this localized-delocalized transformation in terms of critical density nc of fermions (or 
' equivalently, in terms of critical interparticle distance Rc = rifj , where D is the system dimensionality) and the 

■4— ^ ' effective Bohr- radius of the atomic states formed at the transition, in the form n^^^^as ~ 0.25. On the other 
^ ' hand, Hubbard^ and others'* formulated the criterion in the form that at the transition the Coulomb repulsive energy 
^ i among the particles (of magnitude U) is equal to their kinetic (band) energy characterized by the bare bandwidth W , 
' , U :^W. These two criteria are regarded as equivalent. 

^ , The insulator-metal transition is also well defined as a phase transformation in thermodynamic sense. ^ As close 
Q ■ to IMT the band energy (negative) is almost compensated by the repulsive Coulomb interaction, the much smaller 
Q [ thermal or energy can drive the metallic system towards the state with localized spins. The resultant phase diagrams 
contain first-order pressure-temperature transition line with a terminal critical point at temperature Tcr > 0. The 
I [ existence of the critical point was confirmed experimentally quite recently^ and shown to represent properties of the 
^ • Van der Waals liquid-gas point in the D — 3 case of V2 03:Cr, (see Limelette et ali^). The basic question remains 
(~ — ' under what conditions a quantum critical point (QCP) appears at Tcr = between metallic and insulating phases, as 
, its existence marks explicitly the fundamental boundary between the atomic and the condensed (delocalized) quantum 
xj • states. The appearance of such QCP is usually obscured in real systems by the presence of antiferromagnetism on (at 
, least) one side of the transition.^ Nonetheless, the answer to the above question even without inclusion of magnetic 
^sj ' ordering would delineate the basic characteristics of the (quantum) critical point (at T = 0), which is of non-Landau 
characterA 

As IMT involves a drastic change of the nature of states individual carriers, the natural and not addressed so far 
question is how their wave-function changes when approaching the transition from either side. The parameterized- 
model approaches^i^ leave the Wannier functions determining the parameters as fixed. On the other hand, the original 
Mott approach introduces intuitively the concept of an emergent atomic state at the instability of the metallic state. 
^ . Thus the missing question is: how the correlations and the single-particle aspects are interrelated microscopically? 

We provide the explicit description of electronic states close to IMT in which both the interelectronic correlations 
and the single-particle Wannier wave functions {wi{r)} are treated on the same footing. As the result, we introduce 
a singular behavior of the single-particle (Wannier) wave function renormalized by the correlations. This nonanalytic 
behavior (and associated with it unique scaling laws) are only possible when the wave function is calculated a posteriori, 
i.e. in the correlated state, not before the correlations are taken into account. Such an approach is indispensable in 
the situation when the single-particle and interaction energies are comparable or the interaction is even dominant. 
In this sense, our analysis complements that performed within either LDA-|-DMFT^ or LDA-fU^S methods, which 
contain parameters characterizing electronic correlations introduced after the LDA calculations have been carried out. 
We should also underline, that our method was applied earlier to both correlated nano— and macro-^^ystems and 
is free from the double counting of the repulsive Coulomb interaction. Furthermore, the microscopic parameters are 
calculated explicitly within the mutually consistent procedure. Therefore, even though the presented below results 
concern model results, they should be regarded as an essential ingredient to be implemented in modeling real systems 
near IMT. 
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We start from extended Hubbard model with inclusion of intersite Coulomb interaction and ion-ion interaction as 
represented by the Hamiltonian— li^ 

H = el^^ + ^'^^•y "-i" -^U^rii^ n,i (1) 

+ - ^ Kij Sui Suj + Vion-ion , 

where the first term describes the effective atomic energy with 

- + ^Y. {^^J +'2/R,,) , (2) 

the second represents the hopping between the nearest-neighboring sites, the third the Hubbard (intraatomic) and 
the fourth is a part of intersite Coulomb interaction, with Srii = 1 — rii being the deviation from the integer electron 
occupancy rii = 1, and Kij - the intersite Coulomb interaction. The term (2/Rij) in Q expresses the classical 
Coulomb repulsion (in atomic units) for two ions separated by the distance Rij. 

In the Mott-Hubbard insulating state the hole-hole correlations are absent i.e. ((1 — ni){l — rij)) ~ 0. The role of 
these correlations is indeed neglible if the fundamental correlation function (P = {nii^ni{) vanishes for U ^ W . 

The microscopic parameters of this model are expressed via the Wannier functions {z«i(r)} = {w{y — Rj)}liii2 
as follows: ea = {'Wi\Hi\wi) , Uj = {wi\Hi\wj), U = {wiWi\Vi2\wiWi) , and Kij = {wiWj\Vi2\wiWj), where Hi is the 
Hamiltonian for a single particle in the system, and V12 represents interparticle (Coulomb) repulsion. The Wannier 
fimctions are expressed in terms of adjustable Slater functions, i.e. Wi(r) = /34'i(r) — 7X^^=1 ^ji^) where z is the 
number of nearest neighbors, (3 and 7 are mixing coefficients, and ^'i(r) = (a'^/7r)^/^ exp(— a|r — Ri|) is the Is Slater 
function centered on the site i = Rj. In the concrete calculations, they are represented by adjustable Caussians (e.g. 
ST0-7G basis in dimensions D = 1, 2, and 3) and where discussed before i^^i^^ 

As said above, the fundamental principle behind our approach is that the wave functions {wi{r)} are treated on 
the same footing as the diagonalization of ([T]) in the Fock space. Such diagonalization is possible in an exact manner 
for nanoscopic and infinite Hubbard chain {D = 1) systems.'" However, also the Gutzwiller wave-function (GWF) 
and the Gutzwiller-ansatz (GA) approximations lead to close results in the latter case (see below) provided the wave 
functions {wi{r)} are properly readjusted variationally in the correlated state to achieve the ground-state energy as 
a global minimum for a given lattice parameter R. In other words, the energy of the correlated state is readjusted 
iteratively multiple times by adjusting the wave function size a^^ and in effect, also the microscopic parameters. As 
a result of such iterative approach, we obtain the renormalized wave function and the microscopic parameters, as well 
as the ground-state energy, and most importantly, scaling of the physical properties, all as a function of R. Below we 
analyze first the results obtained for GA for three-dimensional cubic lattices, before testing their validity for D = 1 
and 2 situations. The main emphasis is lead on the novel singular scaling properties near the Mott-Hubbard critical 
spacing R — Re- 
in Fig. 1 we plot the exemplary Wannier function centered at R^ = for simple cubic (SC) lattice for R < Rc 
(dot-dashed line) R = Rc (solid line) , and R 3> Rc (dashed line) . One observes immediately that the wave function 
for R > Rc is more extended than that for R = Rc and with increasing R gradually reaches its atomic value oq. To 
study this nontrivial effect in detail, we have plotted in Fig. 2 the relative inverse size 6a/a = \a{R) — a{Rc)\/a{Rc) 
as a function of relative lattice spacing 5R/R = {R — Rc)/Rc in the regime R > Rc- A clear da/dR discontinuity 
(cf. inset) is observed at i? = Rc for all cubic lattices, a rather unique and unexpected feature, which is completely 
absent for the case with bare (unrenormalized) wave functions. Note that the wave function is the narrowest at the 
Mott-Hubbard transition at which « (d^ > in the metallic phase). To determine the universality of the behavior 
we have replotted in Fig. 3 the result of Fig. 2 for i? > Rc in a doubly logarithmic scale. One observes a clear power 
law scaling 6a/a ~ [{R — Rc) / RcY , with s = 0.96 ± 0.1 for not too large R- The results for R < Rc exhibit also a 
similar scaling property for R < Rc, with the exponent s slightly lower (s cri 0.93) and weakly dependent on the type 
of cubic lattice selected. This means that the results in the metallic phase show a lesser degree of universality. Note 
that this power-law scaling describes the singular behavior of the wave-function size as = 1/a. 

Lattice dependence of the scaling with SR/R is also observed for the ground state energy, as shown in Fig. 4. The 
approximate dependence ~ i?* of Ec when approaching the atomic limit can be attributed to the dominant role of 
the Coulomb repulsion between the electrons localized on neighboring ions. 

Both the Figs. 2 and 4 demonstrate the appearance of a simple scaling properties of R^^ type. This dependence 
could be seen in the pressure dependence of the orbital size when close to the critical value R = Rc- However, the 
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FIG. 1: (Color online) Exemplary shapes of the Wannier function centered at the site Ri = for simple cubic (SC) lattice and 
for three lattice spacings marked: R < Rc (dot-dashed line), R = Rc (solid line), and R > Rc (dashed line), ao is the Is Bohr 
radius. 
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FIG. 2: (Color online) Relative inverse wave function size 5a/a as a function of relative lattice parameter SR/R for simple 
cubic (SC), body centered cubic (BCC), and face centered cubic(FCC) lattices. The Mott-Hubbard transition is marked by 
the vertical dashed line. Inset: derivative da/dR vs. 5R/R with a singularity at SR = R~ Rc ~ 0. 



relative changes of a (cf. Fig. 2) are rather subtle and the question is if they can be observable in the present day e.g. 
scanning tunneling observations of atomic orbitals. For that purpose we have plotted in Fig. 5 the maximal value 
wo{0) of the Wannier function relative to that at i? = Rc- The occupancy depletion reaches up to 20% of the peak 
value upon change of 20 — 25% of the lattice constant. Also, one observes the behavior of the same type, as that for a 
in Fig. 2. However, a clear scaling takes place only for R > Rc (see inset), where it assumes the shape ~ (R/Rc)^, 
with u = 0.92±0.01, not much different from that displayed in Fig. 3. On the basis of the results one can say that our 
approach provides evidence for long-range effects in the single-particle wave function which are pronounced already 
in the Mott-Hubbard insulating state. In other words, the Hubbard split subband picture"^ of the Mott insulator is 
concomitant with a strong renormalization of the wave function characteristics near the transition point. 
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FIG. 3: (Color online) Detailed scaling (5a/a vs. 5R/R in the doubly logarithmic scale. The straight lines represent the fitted 
curves ~ {R/RcY . 
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FIG. 4: (Color online) Same scaling as in Fig. 3 but for the ground state energy; the fitted straight-line characteristics are 
listed. 



Finally, we have also examined the role of lattice dimensionality on the singular behavior of the wave function. 
For that purpose we plot in Fig. 6 the inverse wave-function size vs. R/a^ for the linear chain within the modified 
exact Lieb-Wu (LW),— Gutzwiller- wave-function (GWF),— and Gutzwiller-ansatz (GA) ^^i^^'^^ solutions. None of 
the methods provides in £> = 1 case (and GA for D = 2) the type of singular behavior observed in Figs. 2 and 4 for 
D — 3 case. What is more important, all three methods provide quite similar results for D ^ 1. This circumstance 
gives us with some confidence, that the modified above GA method for D = 3 cubic lattice bears physical relevance 
to the problem considered as it removes, among others, a spurious IMT transition for D — I lattice. Therefore, it 
appears that only correlation effect treated in GA combined with the singular behavior of the wave function provide 
a correct IMT characterization. Furthermore, in the inset to Fig. 6 we plot the Wannier function maximum vs. R/uq 
for linear chain (CH), square (SQ), and triangular lattices (TR), respectively. The results for _D = 2 are obtained for 
GA solution. Again, no critical behavior is observed. On the basis of these results we draw the conclusion that within 
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FIG. 5: (Color online) Same as in Fig. 2, but for the maximum wo(0,0,0) of the Wannier function. Inset: Detailed scaling of 
the Wannier function maximum (of same type as in Fig. 3). Fitted straight line represents the function ~ {R/ Rc)^'^^ ■ 




FIG. 6: (Color online) Inverse size of the renormalized wave function for linear chain (CH) vs. R/ao within the exact Lieb-Wu 
(LW), the Gutzwiller- wave-function (GWF), and the Gutzwiller-ansatz (GA) solutions. No singular behavior is observed at 
any R. Inset: Wannier-function maximum ti;o(0) vs. R/ao for linear chain (CH), square (SQ), and triangular lattices (TR), 
with no discontinuity detected in either of the cases. 



our approach the critical behavior appears only for D — i. The striking coincidence of these results to the critical 
dimensionality {D > 2) for the onset of Landau-Fermi liquid stability should therefore be mentioned.— 

One additional basic feature of our approach should be mentioned. Namely, in Table I we list microscopic charac- 
teristic of the Mott-Hubbard transition. These results intercorrelate nicely with the Mott criterion in the following 
manner. The carrier concentration is nc — ^/Rc ^'^^ ^/^c ^'^^ BCC, and for FCC. Thus nc = u/Rq, with 

n — 1,2, and 4, respectively. Therefore, the Mott criterion takes the form: nj/"^ • as — n^^'^ /{RcctC') — 0.20, 0.26, 
and 0.33, for SC, BCC, and FCC lattices, very close values to that obtained originally by Mott.— This connection 
provides additional argument of our approach. 
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TABLE I: Microscopic parameters at the critical point for GA solutions for the lattices specified. 



Struct. 


a/ao Rc/o,Q 


{U/W)c 


SC 


1.099 4.236 


1.337 


BCC 


1.109 4.384 


1.080 


FCC 


1.128 4.351 


0.880 



In summary, we have extended the standard treatment of the Mott-Hubbard transition by incorporating into the 
consistent scheme the wave function renormalization and demonstrating its singular behavior. On the basis of our 
results one can expect also important effect of temperature and magnetic ordering onto this behavior. The long-range 
effects on the wave function shape develop already in the insulating phase and become the strongest at the transition 
by narrowing the atomic wave-function down by about 10%. The wave-function characteristic length plays the 
role of a coherence length for a single fermion in the sea of all other particles. Such an approach can also be extended 
to orbitally degenerate 3(i along the lines proposed earlier.— 
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